Stretching of a single-stranded DNA: Evidence for structural transition 
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Abstract 

Recent experiments have shown that the force-extension (F-x) curve for single-stranded DNA (ssDNA) 
consisting only of adenine [poly(dA)] is significantly different from thymine [poly(dT)]. Here, we show that 
the base stacking interaction is not sufficient to describe the F-x curves as seen in the experiments. A 
reduction in the reaction co-ordinate arising from the formation of helix at low forces and an increase in 
the distance between consecutive phosphates of unstacked bases in the stretched state at high force in the 
proposed model, qualitatively reproduces the experimentally observed features. The multi-step plateau in 
the F-x curve is a signature of structural change in ssDNA. 

PACS numbers: 64.90.+b,36.20.Ey,82.35.Jk,87.14.Gg 
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FIG. 1: (Color online) The schematic representation of ssDNA. Two types of stacking interaction may arise 
namely inter strand shown by the solid line and intra-strand shown by dashed line. In Fig.b we also show 
the steric repulsion among the adjacent bases. The small black circle indicates that one end of the ssDNA 
is kept fixed while a force may be applied at the other end (light blue circle). 

The last ten years have witnessed a revolution in Single Molecules Force Spectroscopy (SMFS) 
experiments involving the manipulation of single biomolecules. These experiments allowed de- 
tecting inter and intra molecular forces and its influence on the properties observable at the chain 
length of biomolecules . New challenges have emerged when the semi-microscopic changes 

in the monomer (nucleotide) found to influence the elastic property of the single-stranded DNA 

nn 

(ssDNA). Attempts have been made to monitor the force-extension (F-x) curve [5|, l6[ of RNA 
and ssDNA consisting of only one type of nucleotide. It was found that the elastic properties of 
ssDNA made up with adenine [poly(dA)] are significantly different from thymine [poly(dT)] (or 
uracil [poly(rU)]). It is known that base stacking is strongest among adenine and weakest among 
thymine and uracil Qj. As a result, base stacking favors a parallel orientation of consecutive bases 
in poly(dA), but not in poly(dT) or poly(rU) [a,|2(]. Notably, theoretical models developed for such 



studies do not include orientation of bases in their description 
limited picture of the system. 
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1 131 ] and hence provide a 



The F-x curves for poly(dT) and poly(rU) show the effect of entropic elasticity, whereas poly(rA) 
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5;| which were found to be absent in earlier experiments [l , 
In a recent experiment, Ke et al. [a] found the existence of multi step plateau in case of poly(dA). 
The first plateau occurs at a force 23 ± 1 pN and overstretched the nucleotide by ~ 74% which has 
been predicted by the model proposed by Buhot and Halperin \)\ . This prediction has also been 
observed in case of homopolymeric RNA [a]. A plateau obtained in case of ssRNA (poly(rA)) has 
been explained on the basis of seven parameters 5[ in the Zimm-Bragg model 



171 ] . The qualitative 



understanding of the first plateau has been achieved in terms of the unwinding of helical structure 
of poly(A) arising due to base stacking. This transition appears to be weakly cooperative, but 



181 ] . Moreover, the 



needs further attention from the statistical mechanics point of view 
second plateau which occurs at a force 113 ± 1 pN and over stretches ssDNA by an additional 
16% Jf! was not predicted by the model proposed by Buhot and Halperin [ltl . Though Ke et al 
conjectured that the second plateau is associated with structural transition [191 ] . but high level ab 



0,0. 



initio quantum mechanical calculation did not support it 

We adopt a more realistic model of DNA where directionality of bases has been included apart 
from excluded volume and non-native base pairing interaction to study the F-x curve. In order to 
study the consequence of orientation of bases, we allow bases to orient along phosphate bonds in 
the model Q|. If two bases are parallel to each other, we say bases are stacked and associate an 
attractive interaction e between them. The issue of cooperativity can also be studied exactly in this 
model by associating a cooperativity factor (<r) between stacked and unstacked domain In this 
paper, we will show that only base stacking and cooperativity factor are not sufficient to describe 
the F-x curve, but further semi-microscopic modifications are needed. Our results based on exact 
enumeration provide unequivocal support that the multi-step plateau is the result of structural 
changes in ssDNA, a field which warrants further studies. 



The model discussed above 21| is general and can be defined in any dimension, though for 
computational limitations, we restrict ourselves in two dimension. We consider a linear polymer 
chain consisting of either (A) or (T) which are described by self-avoiding walks (SAWs) on a square 
lattice. The bases are associated with the link between two monomers as depicted in Fig. 1. 

The present model is a coarse grained one where the microscopic distinction between nucleotide 
A and C has been ignored. Moreover, Seol et al. [f| have shown that the normalized F-x curve for 
poly (r A) and poly(rC) have similar elastic response. In order to study the multi-step plateau [fj, 
we consider following cases: (i) In absence of the stacking interaction, bases are free to orient and 
the ssDNA (or RNA) behaves like a flexible chain. This will correspond to the case of poly(dT) 
[or poly(rU)], (ii) due to the stacking interaction, bases can stack at low temperature which may 
represent the case of poly(dA) [or poly(rA)] With these constraints, we enumerate all possible 
conformations of the ssDNA. The partitions function (Zjsi) of the system can be written as 



Z N = C(N 1 ,N 2 ,x)e^ Nie e l3Fx e- pN2 ^, (1) 

all walks 

where j3 = -^f- Here, e w , T, F and ks are the wall energy, temperature, applied force and the 
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FIG. 2: (a) F-x curve for the model ssDNA for different values of e. e = corresponds to poly(T) (or 
poly(U)), which shows entropic response e ^ represents the case where stacking play a role in the 

formation of helix {e.g. adenine); (b) Experimental F-x curve [fl for the poly(rA). Here F is in pN. (c) 
Figure shows the F-x curve obtained from the model which includes formation of helix (dotted line). In Fig 
a and c force and extension arc in the reduced unit. 



Boltzmann constant respectively. C(Ni,N2,x) is the number of distinct conformations of walk of 
length N. N\ and N2 are the number of intra stacked bases and the number of walls between 
stacked and unstacked region of the chain whose ends are at a distance x apart. a(= e~^ Ew ) is the 
Boltzmann weight for the wall energy which is termed as a cooperative factor. The value of a in 
between 0.5 to 1 corresponds to weakly cooperative, while a — > referred as strongly cooperative 
17]. Time required to enumerate these conformations increases as ji N 2 N because of extra degree of 
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FIG. 3: (a) The Schematic representation of poly(T). (b) In case of poly(A), bases can stack and form a helix 
like structure. Inset shows the base stacking for specific segment; (c) shows the unwinding of helix where 
consecutive phosphate distance remain same; (d) Here the chain is in the stretched state where few bases 
are stacked. Flipping of bases in presence of high force shows the extension in the consecutive phosphate 
i.e a' > a. 
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221 ] . Here (i is the connectivity constant 



freedom associated with bases (2 possible orientation) 
of the lattice. The model also takes care of steric effect among adjacent bases and conformations 
like the one shown in Fig. 1(b) have not been incorporated in the partition function. The quantity 
of experimental interest i.e the reaction co-ordinate in this case is the extension in the chain, 

which can be obtained from the expression (x) = -^J2ni,N2,x x C(Ni, N2,x)e l3Nie e l3Fx e~ l3N2ew . 

In the following, we set fee = 1 and plot the F-x curves for different values of e in Fig. 2(a). 
In absence of stacking interaction (e = 0), we reproduce the entropic response of the chain as seen 

n q 

in case of poly(rll) (or poly(dT)) [5|, laj. In Fig. 2(a) one can see that the curve for poly (r A) 
shifts to the right of poly(rU). Moreover to keep extension say at 0.8, one requires less force for 
polyjrA) compared to poly(rU). In Fig. 2(b), we plot the experimental curve obtained by Seol et 
al. [5(. It is evident that curve for poly(rA) shifts to the left of poly(rU) and to keep extension 
say 1200 nm, one requires almost the double force than the poly(rU). This indicates that apart 
from the stacking interaction and cooperativity there must be some other mechanism involved in 
the elasticity of poly (r A). It was proposed [l(| that since base stacking favors formation of helices, 
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FIG. 4: (Color online) (a) Experimental F-x curve Q for the poly(dT) [dotted] and poly(dA) [dashed] of 
ssDNA. Solid line is for poly(dA) using mFJC (b) Figure shows the F-x curve where the formation of 
helix as well as increase in consecutive phosphate distance between bases has been taken into consideration 
(force and extension are in the reduced unit). One can see that inclusion of such considerations qualitatively 
reproduce the F-x curve obtained by Ke et al (| and exhibits the multi-step plateau; (c) Figure shows the 
probability distribution curves for poly(dA). 



the effective end-to-end distance (R') decreases as shown in Fig. 3(b). Since, we have the exact 
information about the number of stacked bases (N±) participating in the formation of helix, the 
reduction in len gth can be calculated exactly in the proposed model. Following the procedure 



adopted in Ref. 



16[|, we reduce the reaction coordinate x by a certain factor proportional to Ni 



i.e. x' = x — Nia. The value of proportionality constant a = b/a « 0.63 has be obtained by setting 



5.9A° and b = 3.7 A° [5|, LLg]. Here a and 6 are inter phosphate distance and rise in length 



231 ] . has been 



of helix per nucleotide respectively. The modified F-x curve, with this constraint 
shown in Fig. 2(c). Surprisingly this consideration not only shifts the curve to the left as seen in 
expert (earlier shifty to tl . right in absence of helix), hut also q ualitativel y reproduces the 
F-x curve obtained by Seol et al [5J . 

In Fig. 2(a), we show that the force-extension curve obtained by using cooperativity parameter 
a = 0.51 and e = 0.6, which overlaps with a = 1.0 (no cooperativity) and e = 0.8. Therefore, 
stacking interaction reproduces the F-x curve for the low cooperativity. In fact, we find that the 
nature of curve remains almost same upto a = 0.71. This is in accordance with earlier studies 
where it was also found that stacking decreases the cooperativity of melting of homopolymeric 
DNA [24 1 



When cooperativity increases, flatness of the curve [Fig. 2(c)] also increases. The reduction 
in reaction coordinate (x' = x — N\a) along with a modifies Eq. (1) which gives rise such effect. 
From Fig. 2(c), we also note that the curve with cooperativity (a = 0.51) crosses the curve without 
cooperativity (<r = 1) at the extension (~ 0.7). We identified this as a crossover point where chain 



goes from the extended to the stretched state 



21. 
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requires 



Existence of the second plateau [Fig. 4(a)] at higher forces as seen in case of poly(dA) 
further refinement in the model. At such a high force (~ 113 pN) chain must be in the stretched 



state 



221 ] as shown in Fig. 3(d). It was shown 
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26 1 that there may be increase in distance 

□ 

(~ 18%) between consecutive phosphates which drive backbone to a new torsional state [191 ] . 
However 1 — 2 kcal/mol energy difference between two states could not reveal which conformation 
(c3'-endo or c2'-endo) poly(dA) will take under equilibrium condition. In order to see whether 
the second plateau is a result of structural transition, we calculate the number of unstacked bases 
exactly in our model which may participate in this transition. We increase the reaction coordinate 
by a factor proportional to the number of unstacked bases. The proportionality constant for this 
case has been obtained from Ref. [ijj]. The modified F-x curve depicted in Fig. 4(b) shows 
the multi step plateau which is qualitatively similar to the experimental one f|. This implies 
that further increase in the extension (second plateau) is associated with structural transition. 
Moreover, curves for different values of a > 0.51 almost overlaps with the curve for a = 1 and 
hence the transition appears to be weakly cooperative. 

To get the enhanced understanding of the structural changes, we study the probabil- 
ity distribution [P(N±)] of poly(dA) from the following expression (for a = 1): P(N±) = 
~Zj;J2N 1 ,xC'N(Ni,x)e l3Nie e l3Fx . In Fig. 4c, we have shown the probability distribution for stacked 



bases P(Aq) for different values of force at a fixed temperature T = 0.3 for the poly(dA). The 
maxima of P(iVi) for F = occurs at Ni ~ 15. This represents the situation where most of the 
bases are stacked (helix) and domains are randomly oriented. Slight increase in the force aligns 
the helix along the force direction and thus the number of stacked bases increases slightly. This 
can be seen from Fig. 4c where maxima of P(Ni) shifts toward the right. For low forces, thermal 
fluctuations are too weak to unstack the bases in the strand. We observe that there is emergence 
of a new peak around the value N± = 8 and decrease in the peak heights around Aq = 15 at the 
intermediate forces (F = 0.2 to 0.4). The first plateau took place at F = 0.365 as shown in Fig. 
4(b). At this force the height of the both peaks are found to be equal. Two peaks of equal height 
in the probability distribution curve shows the signature of co-existence of two phases (helix and 
coil). This is analogous to liquid-gas transition. Here force represents the pressure and extension 



is analogous to volume [18J. In this region the number of stacked bases decreases which reflects 
the unwinding of helix. This gives relative extension (~ 70%) in the backbone. Further rise in 
force (above F = 0.5) shifts the maxima of distribution curve toward left. This reflects that the 
number of stacked bases in the stretched state decreases gradually and gives additional increase in 
extension which is about 16% as shown in Fig. 4(b). 

In order to resolve the issue associated with structural transition in poly(dA), it is essential 
to calculate the change in energy associated with this transition. Ke et al obtained the total 
energy (3.6 ± 0.2 kcal/mol per base) from the area under F-x curve for adenine (Fig. 2(b) of 
Ref. a]) and interpreted it as base stacking energy. In fact total area under the curve [poly(dA)] 
should be attributed to entropic, enthalpic, structural and elastic contribution. However, in case 
of F-x curve of thymine [poly(dT)], the contributions come from entropic and elastic part only. 
Therefore, at low force (neglecting the elastic contribution) one can get the rough estimate of 
enthalpic contribution (0.6 kcal/mol per base) associated with coil-helix transition corresponding 
to first plateau by subtracting area under II from I below extension 0.58. This value is in agreement 
with the known values available in the literature However, at high force, elastic contribution for 
thymine is different than the adenine and hence one cannot get the energy required for structural 
transition from the same curves. This may be visualized that area under thymine is more than 
the adenine. By using the same persistence length for adenine [la ], one can get F-x curve from 
modified freely jointed chain (mFJC) which contains only entropic and elastic contributions which 
is shown in Fig. 4. The energy (1.2 ± 0.2 kcal/mol per base) associated with structural transition 
above the extension 0.58 has been obtained by subtracting the mFJC curve from the experimental 



F-x curve. This value is in excellent agreement with quantum calculation [2df | . 



We have studied effect of cooperativity in the frame work of exact enumeration. Our analysis 
exactly showed that stacking reduces the cooperativity [Fig. 2(a)] and F-x curve for a upto 0.26 
can be reproduced by stacking alone. For high cooperativity (a < 0.26), stacking is not sufficient 
enough to compensate cooperativity factor. However, formation of helix reduced the reaction 
coordinate which gives an additional factor (aN\) in the partition function. We find that this term 
has significant impact on cooperativity which makes plateau broader in case of adenine. Response 
of cooperativity on the applied force below and above the crossover point is also evident from Fig. 
2(c). 

The exact results on short chain of a new semi-microscopic model presented here provides 
unequivocal support for the structural transition. This has been substantiated by the correct 
analysis of experimental curve along with the mFJC curve that gives the required energy associated 
with structural transition which is in good agreement with quantum calculation. This also resolves 
that c3'-endo puker conformation is most likely in the equilibrium condition and at high force 
poly(dA) goes from c3'-endo to c2'-endo conformation. 
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